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the slow-roll parameters. In the models where the inflaton field is a charged scalar field 
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oscillations of the gauge field and its coupling to inflaton can have interesting effects on 
preheating and non-Gaussianities. 
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I. INTRODUCTION 

Cosmic Inflation proved to be a successful theory of early universe and the mechanism of struc- 
ture formation The simplest modes of inflation are based on a scalar field coupled minimally 
to gravity. The scalar field potential is flat enough to allow for about 60 e-foldings or so to solve 
the flatness and the horizon problem of the standard big-bang cosmology. Precision data can dis- 
tinguish amongst different inflationary scenarios based on their predictions for the spectral index 
of curvature perturbation, the amplitude of gravitational wave or the degrees of non-Gaussianities 
induced on cosmological perturbations. 

There have been considerable interests on primordial anisotropies both observationally and 
theoretically. Observationally, there may be some indications of the the statistical anisotropy of 
the comic microwave background (CMB) [2] although the statistical significances of these findings 
are under debate . On the theoretical sides there have been many attempts to construct 

models of inflation with vector fields or gauge fields which can create sizable amount of anisotropy 
on curvature perturbations These mechanisms can pro vide a seed of anisotropies at the 



order of few percent which may be detectable on CMB [171 ]. [la], 191421], 

There have been different approaches to implements vector field in models of inflation where the 
vector field breaks the gauge symmetry explicitly. One fundamental problem in these models, as 
demonstrated in [22], is the appearance of ghost which render the system unstable and physically 
unacceptable. Therefore, it is crucial that the vector field is protected by a gauge symmetry so the 
longitudinal mode of the vector field excitations is not physical. On the other hand, because of 
the conformal invariance of models with gauge fields, any excitation of gauge field during inflation 
is diluted and can not seed the desired anisotropies. Therefor it is essential that one breaks the 
conformal invariance while keeping the gauge symmetry explicit. This approach was employed in 
different contexts in 



3-30 



In this paper we would like to study different inflationary models where there is a non-zero 
background U(l) gauge field, A^, coupled to a complex scalar field. The charged scalar field can 
be either the inflaton field or the waterfall field of hybrid inflation. Furthermore, in order to break 
the conformal invariance and produce large enough anisotropies, as explained above, we assume 
that the gauge field has a time-dependent gauge kinetic coupling with the kinetic energy in the 
form of =£j&F l j V Fi a '. To be specific, we consider three following models: 

1. Symmetry breaking hilltop inflation: In this model inflaton is a charged complex scalar field 
with the symmetry breaking (Mexican hat) potential V = j (\4>\ 2 ~ ^x~l ■ ^ e wor ^ with 
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the gauge kinetic coupling f(cft) = yj^j with the appropriate parameter p. In the presence 
of a U(l) gauge field, this symmetry breaking potential is theoretically well motivated. 



2. Charged hybrid inflation: This is the standard hybrid inflation [3lf] where the inflaton field 
(ft is real but now the waterfall field tft is charged under U(l) gauge field. The potential 
is v = \ (|Vf - ^-Y + ^(ft 2 \ift\ 2 + ^(ft 2 and the gauge kinetic coupling is f{<ft) = ' " 



where C = M/ g is the critical value of inflaton field at the time of waterfall phase transition. 
3. Chaotic inflation with the potential V = ^\(ft\ 2 where the inflaton field (ft is charged under 



the U(l) field. The case with no gauge coupling (real (ft) was studied in 24|. As in 24| we 
work with the gauge kinetic coupling f(4>) = e c ^ 2 / 2M p with the appropriate coupling c. We 
generalize their results to our case with a non-zero gauge coupling. 

The rest of paper is organized as follows. In section [II] we study our set up where the inflaton 
field is charged under the £7(1) gauge field with a time depending gauge kinetic coupling. 
In section [III] we study the symmetry breaking potential in details and calculate the level of 
anisotropy during inflation. In section IIVI we concentrate on hybrid inflation. We will see 
that the presence of the gauge field plays important role at the end of inflation and dynamics 
of waterfall phase transition. In section [V] we briefly study the case of chaotic inflation where 
now the inflaton field is charged under the U{\) field and compare our results with those of 
. Conclusion and brief discussions are given in section IVII 

II. BACKGROUND EQUATIONS 

Here we present the action for the cases of symmetry breaking potential where the inflaton field 
(ft is charged under the U(l) gauge field with a ^-dependent gauge kinetic coupling f 2 ((ft). The 
action and the background equations for the other cases can be obtained accordingly. 

As in 1 261 the action is 



/ 



S = d A Xy/—g 



2 1 t2t.... 

(1) 



M z I f 
-fR - -D^D^ - LJZL FiafF v» - V ((ft, 



where Mp 2 = 8vrG, for G being the Newton constant and the overline represents the complex 
conjugation. The covariant derivative is given by 



Dfj,(ft = d^(ft + ie(ftA fl 



(2) 
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where e is the dimensionless gauge coupling of to <fi. As usual, the gauge field strength is given 
by 



Fav — V nA v — S/ u An — dnA v — d v Au . 



(3) 



We work with potentials which have axial symmetry where V and /(</>) are only functions of 
4>(f> = \4>\ 2 . It is more instructive to decompose the inflaton field into the radial and angular parts 



<f>{x) = p{x) e i9 ^ 



(4) 



so V = V(p) and f 2 (<f)) = f 2 (p)- As usual, the action ([!]) is invariant under local gauge transfor- 
mation 



e 



+ e(x) 



With this decomposition, the action (pQ) is transformed into 
r2 r „2 



S 



d x^J—g 



M' z i n 2 f 2 (n) 



The corresponding Klein-Gordon equations of motion are 



accompanied by with the Maxwell's equation 



d,{V=g~f\p)F^)=er, 



where the current J v is defined by 



J u = p 2 V=g{d u 6 + eA») 



(5) 



(6) 



(7) 
(8) 



(9) 



(10) 



The conservation of from Eq. ([7]) is a manifestation of the axial symmetry imposed on V. 
Interestingly, Eq. ([7]) is not independent from Maxwell's equation, where F^ v being anti-symmetric 
leads to d^d v F^ = = 0. 

Finally, the stress energy momentum tensor, T a p, for the Einstein equation, G a p = 8TrGT a p, is 



la/3 



'-QafiF^F^ + f(p)F a ^ + d a pd pP + 



9a/3 



.(11) 



4 - ' ' ■ ' p*\g\ 

We are interested in the effects of a non-zero background gauge field on the evolution of system. 
To fix the gauge we use the Coulomb-radiation gauge Aq = diA 1 = 0. Since our background is only 
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time-dependent, from the constraint J° = one concludes that 6 = at the level of background. 
The inclusion of a non-zero background gauge field breaks the Lorentz invariance explicitly since 
a preferred direction is singled out in the background space-time. We take our background gauge 
field to have the form = (0, A(t), 0, 0). 

As in [24j] our background metric has the following form 

ds 2 = _ dt 2 + e 2a(t) {e -Mt) dx 2 + e 2a(t) {dy 2 + ^ 2)) {u) 

Here a(t) measures the background number of e-foldings, d represents the background isotropic 
Hubble expansion rate while &{t) measure the anisotropic expansion rate. For a universe with 
small anisotropies, we require that \&/a\ <C 1. 

Assuming that the fields p, A, a and a are only function of t the background equations of motion 

are 



d t [f 2 (p)e a+4a A) = -e 2 p 2 e a+ ^A (13) 
p + 3ap + V p + (-f( p )f p (p)A 2 + e 2 P A 2 ) e~ 2M = (14) 



\p 2 + V(p) + [lf(p)A 2 + ±f A 2 ) e' 2 ^ = 3M 2 (a 2 - a 2 ) (15) 
V (j,) + ( \f(p)A 2 + ^A 2 ) e~ 2M = M 2 (a + 3d 2 ) (16) 



(17) 



6" xr ' 3 

'\f\p)A 2 - ^fA 2 ) e-***> = Mp (3da + a) 

n 

In the limit where e = these equations reduce to those of [24]. One can also check that not all 
equations above are independent. For example, Eq. (|15p can be obtained from the remaining four 
equations. 

From Eq. (|14p . the total energy density, £, governing the dynamics of the inflaton field is given 

by 

i T/ _i_ ^-2a+4a ( 1 f 2f„s X2 , ^ A 2 



8 = '- + V + e"™ ( |/ (p)A 2 + -^A^ . (18) 

Since the second and the third terms above come from the gauge field, we refer to them respectively 
as the kinetic energy and potential energy associated with the gauge field. Using Eqs.(15) and (16), 
the equation for acceleration of the universe is given by 



d + d 2 = -2a 2 + ' 



3Mf 3Mp 



V-±f(p)A 2 e- 2 ^ 



(19) 



Interestingly, the term in effective potential proportional to e cancels out in this expression. One 
also observes that for inflation to take place, corresponding to d + d 2 > 0, one requires that the 
background potential V(p) dominates over the contribution from the gauge field kinetic energy. 
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In the following we are interested in configuration where inflation take places with small 
anisotropies such that 

<5=|?|«1. (20) 
a 

Small amount of anisotropies in background inflationary dynamics may be acceptable assuming 
that they do not impose too much anisotropies on CMB temperature power spectrum. In order 
for anisotropies to be small, the contribution of gauge field to the total energy density should be 
small compared to the background potential. To parametrize this, we define the ratios R\ and R2 
via 

„ A 2 f(p) 2 e- 2a n e 2 p 2 A 2 e~ 2a , . 

2V ' R ^ 2V ■ (21) 

For the contribution of the gauge field energy density into the total energy density to be small we 
require R\ , R2 *C 1. In this limit where the anisotropy is smaller than the slow-roll parameters 
(defined below), from Eq. Q17|) we obtains 

5 ~ - Jfe) . (22) 

One of our goal in this work is to see the behavior of R\,2 during inflation. As we shall 
see, during early stage inflation both R% and R2 are very small and inflation is basically driven 



by the background potential V and one can treat the system as isotropic inflation. As in 
sometime during inflation, R\ rises quickly such that its contribution to the Klein-Gordon equation 
governing the scalar field dynamics can not be neglected. This is an attractor mechanism and as 
we shall see below R\ scales like the slow-roll parameters once the system is in the attractor regime. 
Interestingly, the Hubble expansion rate is still predominantly given by the background potential 
V but one should check that the anisotropy given by Eq. (|17p is under control. One new effect in 
our model is that sometime at the end of inflation R2 becomes comparable to R\ and the gauge 
field oscillates very rapidly. Because of the interaction term e 2 p 2 A 2 , the rapid oscillations of the 
gauge field induce rapid changes in inflaton effective mass, violating the slow-roll conditions and 
ending inflation abruptly. Here we would like to study these three distinct inflationary phases in 
some details. For convenience, we refer to these three stages of inflation as phase one, two and 
three, respectively. However, note hat the third inflationary stage is very short compared to other 
two inflationary stages. 

As a measure of slow-roll parameters and phase change we define the dimensionless quantity e 
and 77 given by 

e = ~ » n = -^- ( 23 ) 

or Sap 
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When the slow-roll approximation holds e, r/ <C 1 and inflation ends when e, rj ~ 1. In our 
anisotropic inflationary models, 77 has sudden jumps which represents the onsets of phase changes. 
In Fig. [T]we have plotted 77 for some parameters value which clearly indicates the jumps in 77. 



III. SYMMETRY BREAKING HILLTOP INFLATION 



We start with an almost isotropic configuration with negligible anisotropics such that the gauge 
field contributions in expansion rate Eq. f)15[) and the Klein-Gordon equation (|14p are negligible, 
corresponding to i?i,i?2 <C e,ij. However, we would like to allow the gauge field kinetic energy to 
increase such that 5 increases towards the allowed observational bounds. 

To be specific, we work with the symmetry breaking (Mexican hat) potential which is theoret- 
ically well motivated in a model with an Abelian gauge field 

Here A is a dimensionless coupling and for the later convenience we have defined the dimensionless 
variable p = p/p where p = M/y/X. The potential has global minima at p = p or p = 1. In 
this picture, we assume inflaton starts at the top of the potential and proceeds towards the global 
minima. As is well known cosmic strings are produced at the end of inflation which can have 
interesting observational effects. 



In 



241 ] the authors studied the simple chaotic inflationary potential V = m <j) /2 for a real scalar 



field. As we shall see, their conclusion about the existence of the attractor mechanism during the 
second phase and the behavior of R\ and 5 will also hold in our case. However, the contribution of 
gauge coupling e via R2 opens up new inflationary phase and the dynamics of the system at the 
end of inflation is quite different than what studied in 24 ] . 



During the inflation, which happens mostly in the hilltop regions of the potential, one may 
approximate potential (JMJ) as 

M 4 M 2 9 

v ^lx--»- (25 » 

For this approximation to be valid, one has to satisfy p <C 1 during inflation. 

Since we are interested in small anisotropies, R±,R2,S 1, the background expansion is given 
as in standard slow-roll inflationary models with 

tf^J^-* (26) 

3MJ, 12AM|, V ' 
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FIG. 1: Here we have plotted rj defined in Eq. (J23J, with A = 2.5 x lO" 13 , M = 3 X 1CT 6 M P , p = 50, 
Pi n = Mp/5 and e = 1. The first phase change, at e-folding a ~ 10, happens when R\ becomes comparable 
to e. The second phase change happens very close to the end of inflation (in about one e-folding towards 
the end of inflation) when i?2 also becomes comparable to e. 



Our numerical investigations show that the phase changes happen when there are sudden 
changes in p in a short period. This can be seen in the plot of rj in Fig. Q] where there are 
two sudden jumps. The first jump corresponds to the transition from phase one to phase two 
where the R\ contribution in Eq. (|14p becomes important. The second jump, very close to the 
end of inflation, represents the transition from the second phase to the third phase. This happens 
when the right hand side of Eq. (|13p can not be neglected and it eventually affects the evolution 
of p in Eq. (|14p . As we shall see the third period is very short and inflation ends abruptly once 
the gauge field starts to oscillate during the third phase. 



A. Phases One and Two 



At the early stage of inflation, the contribution of gauge field in total energy density and the 
scalar field equation is completely negligible corresponding to -Ri,i?2 "C e,rj. Inflation proceeds as 
in standard slow-roll hilltop inflation and in order for the slow-roll condition to be satisfied one 
requires that Mp{V p /V) 2 and MpV pp /V both to be much smaller than one. These in turn yields 
p c ^> 1 where 

M 2 



Pc = 2\M? ' (27) 
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The Hubble expansion rate is given by Eq. (|26p while the scalar field equation in the slow-roll 
approximation is 

P^^TP P^Pine 2a/p °, (28) 

where pi n represents the initial value of the inflaton field. Also here and below, prime denotes 
derivative with respect to a, the number of e-foldings. We use the convention such that at the 
start of inflation a = and the total number of e-foldings measured at the end of inflation is 
a = atf — 60 to solve the flatness and the horizon problem. From Eq. f|28|) the number of e- folds 
as a function of p during the first phase is 

a(p) ~ f In (JL). (29) 

* \ Pin J 

So far we have not specified the form of f(p), the time-dependent gauge kinetic coupling. In 
order for the perturbative gauge theory to be under control we demand that the effective gauge 
kinetic coupling <M(a) = f{p)~ l to be small during inflation and approaches unity at the end of 
inflation for some yet unknown dynamical mechanism, that is g^(af) = 1. This indicates that 
f(p) is a decreasing function during inflation. As mentioned before, we would like the gauge 
field contribution to the energy density to be subdominant but big enough to play some roles in 
anisotropy and scalar field equations. To determine the form of f(p) we note that during the 
first two phases, R2 <C R\ so the terms proportional to e in background equations (fT3l) - (fT71) can 
be neglected. From Eq. (fT3j) one obtains A oc f(p)~ 2 e~ a so R\ scales like R\ oc f(p)~ 2 e~ 4a ~ 
f(p)~ 2 P~ 2pc - Therefore, for the critical coupling f c = (p/p) Pc the gauge field kinetic energy remains 
fixed during the first two phases. As we started with negligible R± in phase one, then it remains 
negligible afterwards, i.e. R\ -C e. In order to increase R\ during the second phase we consider 
the gauge kinetic coupling 

f( P )=(^y =p- p , (30) 

with p > p c such that the gauge field kinetic energy becomes important during the second and 
third stages of inflation. 

As explained above, during the first two phases the right hand side of Eq. (|13p can be neglected 
and using Eq. (|30l) one obtains 

A!={^Ye~\ (31) 

where the dimensionless gauge field is defined via A = A/p and £ is a constant of integration. Note 
that we defined the constant of integration in this way so the subsequent analysis becomes more 
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simplified. Physically, £ is measured by the initial value of Ri at the start of inflation, a = 0. 
Plugging Eq. (j3T1) into Eq. (|2TT) during the first two phases one obtains 



R 



and therefore the initial value of R\ is 



(32) 



Ri 



. ~ E£ (to- ) 2p 
in — g ys.Pinj 



(33) 



Since we demand that R\^ <C 1, the Friedmann equation is still given by Eq. (|26p . Combined 
with Eq. ()30p . the inflaton field equation in the slow-roll limit is cast into 



2\/ 



(34) 



3 

As can be seen from Eq. (I32p . (£,p) 2p is very small during the first phase of inflation so one can 
neglect the last term in Eq. (|34p and the solution is given by Eq. (|28p . For this to take place, 
we need to make sure that at the start of inflation, a = 0, the third term in Eq. ()34|) is indeed 
much smaller than the second term. Using the expression of Rn n given in Eq. (I33p this condition 
is transformed into 



Rl in ^ Pin ■ 
P 



(35) 



As inflation proceeds and p increases the last term in Eq. (j28|) catches up with the second 
term and one should take the effect of this term into account. This is exactly when the gauge field 
contribution into the inflaton equation, Eq. (|14p . becomes important as promised. Eq. (|28p can 
be solved with the answer 



P 



2a 

Pine vc 



1 + 



6(p-Pc) 



\2p 



M.p — Pc)a 
1 Pc 



l/2p • 



(36) 



During the first inflationary phase, the second term in the denominator is much smaller than unity 
and the solution to the above equation reduces to our previous result, Eq. (|28p . The transition 
from the first phase to the second phase happens when the two terms in the denominator above 
become comparable. Defining the first phase transition to take place at a = a±, one obtains 

6(p - p c ) 



a i 



Pc 

4(p - p c 

Pc 

4(p - Pc 



In 
In 



P 2 Pc (£p* 
2(P-Pc) 



\2p 



p 2 R 



1 in 



(37) 
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FIG. 2: Here we plot our analytical solution for p(a), Eq. (|5tH) . shown by the red dashed curve, and compare 
it to the full numerical solution denoted by the solid black curve. The agreement between them is very good. 
The left figure corresponds to e — 1 whereas for the right figure e = 10~ 4 . As argued, the time of first 
phase change, which here is at a\ ~ 10, is independent of the value of e and is well approximated by our 
analytical formula Eq. (|3T[) . All other parameters are as in Fig. [T] 



where to get the final answer Eq. f)33f) has been used. This is an interesting result because the 
onset of the first phase transition is controlled by the anisotropy at the start of inflation, Rn n , 
and the parameter p. As can be seen from Eq. (|37p . the smaller is the value of initial anisotropy 
Rim, the longer it takes for the system to enter the second inflationary regime. We have checked 
that Eq. (|37p gives a good estimate of a\ compared to the full numerical results. 

In Fig. [2] we have plotted our analytical solution for p(a) compared to the full numerical 
analysis. As can be seen they are in very good agreement. Since Ri <C R\ in this period, the 
value of a\ is independent of the gauge coupling e as can be seen from the plots. The left figures 
correspond to e = 1 whereas for the right figure, e = 1CP 4 . 

In order for the phase transition to take place during the 60 observable e- folds, we require 
a.\ < 60. This in turn impose the following lower bound on Rn n 



2(p-p c ) 
R\ in > 5 exp 



240(p - p c 



(38) 



Pc 

This is reasonable, because the smaller is the initial value of anisotropy, the longer it takes for the 
gauge field kinetic energy to become significant to affect the dynamics of the inflaton field. For 
values of p c comparable to p Eq. (|38p can easily be satisfied and the first phase transition takes 
place during the physically relevant window of inflation. 

During the second phase, a > a±, the solution ([36]) quickly approaches to its attractor solution 
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FIG. 3: Here we plot the evolution of the gauge field. The left graph represents In A where the red dashed- 
dotted curve is our analytical solution Eq. (J49j) whereas the solid black curve is the full numerical solution. 
The agreement between our analytical solution Eq. (|4TJ)) valid for the second and third phase and the full 
numerical result is good. Also the change of the slope of In A form the first phase to the second phase 
is in good agreement with our other analytical result, Eq. pt)|) . valid for the first two phases. The right 
graph represents A during the last few e-foldings. The oscillatory behavior suggested by Eq. (|4l?)) is clearly 
seen. The start of the third phase, corresponding to the first peak is well approximated by our analytical 
estimation Eq. (|47|) . All parameters here are as in Fig. [T] 



and 



V Vc 



(39) 



If we plug back this equation into the inflaton equation (|34p . we find that the last term in 
behaves as a constant source term with the magnitude 4(p — p c )/pp c turned on at the time of phase 
change. This explains the kink in p behavior seen in Fig. [2j 

It is also instructive to look into the gauge field evolution in this stage. Plugging Eq. ([36]) into 
Eq. (|3ip one can find an analytic expression for A' valid for both phase one and two with the 
asymptotic behavior 



A' 



{AY exp 

6(P-Pc) e 3a 



(4p-p c )a 

Pc 



OL < OL\ 

a > a\ 



(40) 



This indicates that the gauge field increases like e 3a during the second phase whereas it was 



increasing with a slightly higher rate, exp 



(4p— p c )a 

Pc 



, during the first inflationary stage. The 



change in the slope of the evolution of In A clearly can be seen in Fig. [3l 
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Finally we are in a position to find the form of R\ and 5. During the first inflationary stage, 
plugging Eqs. flU} and (J29]) into Eq. ([H]), yields 



2 2 
5 ~ ~ -Ru n exp 



4(p-p c )a 



(41) 



Pc 

As expected, 5 increases exponentially during the first inflationary stages with the initial amplitudes 
set by R\i n . During the second stage it reaches its attractor value. Plugging Eq. ([39]) into Eq. 
(|2ip . during the attractor regime we have 



^ggi^-/c; (42) 
3 3p 2 



This attractor value is fairly independent of the initial conditions. 



In chaotic model studies in 



24j ] it was shown that R\ follows the slow-roll parameters R\ ~ e. 



Here we show this conclusion also holds for our case. To see this, note that in the slow-roll limit 
e~2R 1 + ^ F . Using Eq. (J39]) in Eq. (JM]) one can approximately find that e ~ 2R\ + 4p 2 £> C/ /p 2 . 
Since p < 1 during the second phase one concludes that 

e > 2R X . (43) 

In Fig. U]we have plotted the ratio R\/e and the anisotropy 5. As can be seen, our analytical 
formulae Eqs. (|43h and (142)) are in good agreement with the full numerical results. As explained 
before, the solution during the second phase quickly reaches the attractor regime where R\ ~ e ~ 5 
and the fraction of gauge field energy density to the total energy density in the Friedmann equation 
is at the level of slow-roll parameter. The attractor phase can clearly be seen from the behavior of 
Ri and 5 in Fig. BJ 

As the gauge field A increases exponentially the effective potential for the inflaton field increases 
as e 2 A 2 p 2 e~ 2a oc e 2 p 2 e 4a and the slow-roll condition quickly terminates at the final stage of 
inflation. This is the third stage where R2 becomes comparable to R\. Below we study this phase 
in some details. 



B. Final stage of inflation 

Now we consider the final stage of inflation when the right hand side of Eq. ()13|) can not 
be neglected. Using Eq. (f39j) in Eq. (fT3|) . the equation of motion for the gauge field can be 
approximated to 

A" - 3A' + fie 4a A = (44) 
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FIG. 4: In left figure we plot ln(i?i/e) ( upper solid green curve) and ln(i?2/e) (lower dashed-dotted red 
curve). The phase change takes place at ai ~ 10 followed by the attractor regime denoted by the almost 
horizontal line where Ri oc e. As explained in the text, i?2 is very small compared to R\ until the very end 
of inflation when they become comparable and inflation ends shortly after that. Right: ln<5 is presented. 
The attractor behavior during the second inflationary stage is clear. All parameters here are as in Fig. [TJ 



where the dimensionless parameter (3 is defined via 

36e 2 (p — p c 



P 



A p 2 Pc £, 2p 

The solution to this equation is in the form of Bessel functions 



3a 

e 2 



VP 



ai J 3/4 I lJ-e 2a j + a 2 Y 3/4 I ^-e. 



VP 



2a 



(45) 



(46) 



where a\ and a 2 are the constants of integration. 

From the form of Eq. ()44p it is seen that the third inflationary phase starts when the last term 
in Eq. (|44p is comparable to the second term. This means that VP e2a2 — 1 where a = a 2 is the 
start of the third inflationary stage. This gives 



1 



a 2 



In 



(47) 



36 e 2 (p — p c 

We have checked numerically that this expression gives a very good estimate for a 2 , the onset of 
transition from the second inflationary stage to the third inflationary stage. Shortly after a > a 2 , 
the argument of the Bessel function exponentially increases and the gauge field starts to oscillate. 
This in turn triggers a sharp increase in the slow-roll parameters e and rj and inflation ends abruptly. 
This can be seen in the plot of rj shown in Fig. [TJ 



15 



For the consistency of our setup we require that a\ < a 2 , i.e. the third inflationary stage 
takes place after the second inflationary stages. Comparing Eqs. (j4"Tl) and (i3T|) . and assuming 
(p ~ Pc) I Pc ~ 1, one requires 

e 2 «Ap 4 £ 2p iW (48) 

Because pi n <C 1 and p>l, this condition can easily be met. 

The condition y^e 2 " 2 ~ 1 indicates that during the final stage of inflation the arguments of the 
Bessel functions in Eq. (149p are bigger than unity while during the first two inflationary stage the 
arguments of the Bessel functions are small. Using the small argument limit of the Bessel function, 
J 3 / 4 (x) ~ x 3 / 4 and Y3/4 (x) ~ x~ 3 / 4 for x <C 1 one concludes that the term containing Y 3/ / 4 decays 
quickly as inflation proceeds and the term containing J 3 / 4 survives in Eq. (|49p . More specifically, 
J3/4 (yT? e 2a /2) ~ (v^5/2) 3 / 4 e 3a//2 and comparing this with Eq. (|40p during the second inflationary 
stage one can fix the coefficient a\ to obtain 

Note that this expression works for the second and the third inflationary stages whereas the formula 
Eq. (|40p works for the first two inflationary stages. We have checked that Eq. (|49p is in good 
qualitative agreement with the full numerical analysis. In Fig [3] we have compared Eq. (|49|) with 
the full numerical result and the agreement between them is good. Also in Fig [3] we have plotted 
the behavior of A for the last few e-foldings. The start of the third phase is when the argument of 
the Bessel function in Eq. (|49|) becomes at the order of unity given by Eq. (|47p . This corresponds 
to the first peak in the plot of A in the right figure of Fig 

It is also instructive to compare R 2 with R\ during the final stage of inflation. As explained 
before, R\ and R 2 measure respectively the gauge field kinetic energy and potential energy com- 
pared to the background inflationary potential. Physically, we expect that during the final stage 
of inflation R 2 rises quickly and becomes comparable to R\ and e. In Fig [4] we have compared R 2 
with R\ and 5. Initially R 2 is very small, but during the final stage of inflation R 2 rises quickly and 
becomes comparable to R\. Physically this means that the inflaton mass receives a time-dependent 
contribution of the form e 2 ^4 2 p 2 e~ 2a and the slow-roll conditions are violated soon after the gauge 
field starts to oscillate. This conclusion is supported in both Fig [3] and Fig [TJ 

Finally it is also instructive to look into the behavior of the inflaton field as a function of the 
the strength of the gauge coupling e. As argued before, during the first two inflationary stages e 
does not play important roles and the evolution of the inflation proceeds as in e = 0. In particular, 
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the position of the first kink, a±, is quite insensitive to the value of e. On the other hand, e 
controls the end of the second inflationary phase, ct2, but only logarithmically. In the right plot 
of Fig. [2] we have changed e by four orders of magnitudes. Correspondingly, 02, and the total 
number of e-foldings changed by 4. This is consistent with our analytical formula Eq. (|47p which 
for e —7- 10 _4 e predicts an increase of e-foldings of (In 10 8 )/4 ~ 4.6. 



IV. CHARGED HYBRID INFLATION 



In the previous example the inflaton field was a complex field and was responsible for the 
symmetry breaking. Now we consider the case where the inflaton field is real, while the symmetry 
breaking is controlled by another complex scalar field, the waterfall field. The action is 



d x\/—g 



M 



1 



(50) 



As explained above, cj> is the inflaton field while ip is the complex waterfall field. 
The potential is as in standard hybrid inflation j3l| 



A 



M 2X ~ 2 



A 



+ 



9 A*\ 



+ 



m 



(51) 



We are interested in the configuration where the potential is axially symmetric and V(ip, if), eft) = 
V(x,4>) where ip(x) = x(x) e* 6 ^). Following the same metric ansatz as in Eq. (|12p and taking 



= (0, A(t), 0, 0) the equations of motion are 



$ + 3 d</> + 0(m 2 + g\ 2 ) - /(<A)/40)i 2 e- 2Q+4CT 
X + 3«X + (\{X 2 - ^) + 9 2 A X + e 2 xA 2 e" 2a+4ff 



\tf + \x 2 + Vfax) + { \f\<t>)A 2 + ^A 2 ) e~ 2 ^ 



2 



2 2 

2 , e X ,2 \ -2a+4o 



r (( , xM -(^/ 2 (0)^ + ^-^) ( 



-e 2 X 2 e a+4a A 





3M| (a 2 - a 2 ) 
Ml (a + 3d 2 ) 
Ml (3d<r + ir) . 



(52) 
(53) 
(54) 

(55) 

(56) 

(57) 



From Eq. (|55p . the total energy density determining the expansion rate of the universe is given 



by 



£ = V(<t>, X ) + e- 2a+4 ° ( \f\<P)A 2 + ^A 2 



(58) 
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The interesting new effect is that the gauge coupling e induces a new time-dependent mass 
term for the waterfall field in the form e 2 e~ 2a A 2 \ 2 ■ This can be seen both in total energy density 
and also in equation governing the dynamics of the waterfall field, Eq. (154j) . As in standard 
hybrid inflation we work in the vacuum dominated regime where the waterfall field is very heavy 
during inflation so % quickly settles down to its instantaneous minimum \ = during inflation. 
In standard hybrid inflation models, inflation ends when inflaton field reaches a critical value, 
<fi = cj) c = where the waterfall field becomes tachyonic and rolls down very quickly to its global 
minimum ip = fi = Mj y/X, = ending inflation very efficiently. In our model to find the moment 
when the waterfall field becomes tachyonic, let us calculate the x field effective mass 

0lx=o=5 2 (^ 2 -^) + eW. (59) 

In the absence of the gauge field, the onset of waterfall field instability is when <p = <j) c . However, 
in the presence of the gauge field the time when the tachyonic instability is triggered is modified. 
Indeed, if either of e or the background gauge field A are very large, then the onset of tachyonic 
instability can be significantly altered and inflation will end before 4> reaches 4> c . This can have 
profound effects on the dynamics of waterfall phase transition and symmetry breaking 
Furthermore, one needs to revisit the question of tachyonic preheating in this case. 
The condition of waterfall phase transition, Eq. (I59p . can be rewritten as 



32, 



3J. 



e 2 - 

cf) 2 + -^A 2 e~ 2a -1 = (60) 

r 

where we have defined the dimensionless fields <f> = (fi/4> c and A = A/(f> c . In this notation, in the 
absence of gauge field the waterfall phase transition happens at <j) = 1. If the gauge field is expected 
to play important roles in determining the dynamic of waterfall phase transition then one requires 
the second term in Eq. (|60p to become at the order of unity at the time of transition. Below we 
will study under what conditions on model parameters this condition can be met. 

As mentioned above, we assume the waterfall field is very heavy during inflation and the po- 
tential driving inflation is 

M 4 1 o 9 

y ^4A + 2 m0 ' (61) 

In order for the inflaton field to be light during inflation so the slow-roll conditions are met we 
need p c 3> 1 where now p c is defined via 



18 



Furthermore, the assumption that the waterfall field is very heavy during inflation requires 
XM 2 p jM 2 > 1. Finally, the condition of vacuum domination during inflation is met if X/g 2 <C 
M 2 /m 2 . 

As in our previous symmetry breaking example we assume the anisotropics are negligible corre- 
sponding to 5 < e throughout inflation so the background expansion is still given by Eq. (I26D . As 
in symmetry breaking example, inflation starts with the isotropic limit where R\ <C e at the early 
stage of inflation. As inflation proceeds, R\ rises quickly and we enter the second phase of inflation 
where the gauge field dynamics affect the evolution of the inflaton field in Eq. (|53p . However, 
unlike the previous example, the final stage of inflation will be very different where now inflation 
ends violently once the waterfall field becomes tachyonic. 

During the first phase of inflation, the inflaton dynamics is 

0' + 2 ± = o _> ~ 4> in e- 2 ^ , (63) 

Pc 

where (pi n is the initial value of the inflaton field. The number of e-foldings is 

— HI)- 

Now we need to determine the form of the gauge kinetic coupling, /(</>), such that R\ rises 
quickly during the second phase of inflation. Since during inflation \ = 0, then the gauge field 
equation ([52]) can easily be solved with A' ~ e~ a f((j))~ 2 . Consequently, R\ scales like Ri ~ 
/ '(4>)~ 2 e~ 4a ~ f (4>)~ 2 (f) 2pc . Therefore, for the critical coupling f c ~ <£P C , R\ remains fixed and the 
fraction of the gauge field kinetic energy to the background energy density remains fixed. As in 
previous example, we allow for the following coupling 

m -(£)'-*, (66) 

with p > p c so the energy density of the gauge field increases as inflation proceeds. As in the 
symmetry breaking case, we expect to enter the attractor regime where R\ ~ e till inflation ends 
via tachyonic instability. With gauge kinetic coupling given by Eq. (i65|) . the gauge field equation 
can be solved easily and 

A' = e- a (i^)~ P , (66) 

where £ is a constant of integration. Plugging the gauge field solution into the inflaton field equation 
yields 

(^y + 0_ _ M^l^yp e -4a = (67) 
Pc 3 
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In comparison to symmetry breaking analysis, it is interesting to note that one can reproduce the 
previous results with the replacements p — > —p and p c — > —p c - Like in symmetry breaking example, 
this equation can be solved easily with the solution during the first stage of inflation given by Eq. 
([63]) while during the second stage of inflation, for a > a\, one obtains the attractor solution 

Furthermore, plugging Eq. ([68]) into gauge field equation ([66]) yields 



A ^ 2( V c) ^> . (69) 

P Z Pc 

As expected, the gauge field increases exponentially. 

We can also calculate the level of anisotropy in this scenario. For the first inflationary stage, 5 
is given as in Eq. (]4ip whereas during the second stage it reaches the attractor value 

^3 El " Sp 2 #M? ' (70) 
As expected, R\ reaches the scaling solution during the second inflationary stage. Like in symmetry 
breaking case one obtains e > 2R\. 

Now we have all the tools to answer our original question that under what conditions the gauge 
field can play a role in triggering the water field phase transition. As explained below Eq. (|60f) . this 
can happen when the combination (e 2 / g 2 )A 2 e~ 2a f is comparable to unity at the end of inflation 
when a = atf — 60 and <j) = 1. Using Eqs. (|68p and (|69p . and noting that by definition = 1 when 
a = ctf, one obtains the interesting results 

S± e \ctf — o22 2 2" V IS -> 

This indicates that for for e <C pg the gauge field does not play important role in triggering the 
waterfall field tachyonic instability and inflation ends as in standard hybrid inflation. However, for 
e 3> pg, then the gauge field shuts off inflation before <p = <j) c and the dynamics of the waterfall 
phase transition, symmetry breaking and tachyonic preheating would be drastically different than 
what happens in standard hybrid inflation. Because of the inhomogeneous end of inflation large 
non-Gaussianities can be produced at the end of inflation in the light of [l^ |. We would like to 
come back to this question in a future publication. 



V. CHAOTIC INFLATION 



As our final example, here we briefly study the case of chaotic inflation. Many of our previous 
results also apply here. For this purpose, we will be brief here only emphasizing our main results 
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and compare them to the results of 24| where they assumed e = 0. To be specific, we concentrate 
on the quadratic potential 



m 



V = —I 



(72) 



The equation of motions are the same as in symmetry breaking example. During the first infla- 
tionary stage, the relevant equations are 



, m 2 p 2 

a ~ — 

~ 6M 2 



3ap + m p ~ . 



(73) 



where p = \cp\. Note that in our slow-roll limit the Friedmann equation above also holds throughout 
the inflationary period. These equations can easily be solved to give 



p 2 ~ pj n - AMpOt. . 



(74) 



Now we determine the form of the desired gauge kinetic coupling. As in previous examples 
we start with an almost isotropic configuration with negligible anisotropics such that the gauge 
field contributions in expansion rate Eq. (|15p and the Klein-Gordon equation (|14p are negligible. 
However, we would like to allow the gauge field kinetic energy to increase such that 5 is within the 
observational bounds. As can also be seen from the solution of Eq. (|13p . in order for the gauge 
field energy density to remain constant during the first inflationary stage the gauge kinetic coupling 
should have the critical form f c ~ e 2a . For the power law inflationary potentials with V oc p n the 
critical gauge kinetic coupling is given by f c ~ e P 2 / nM p_ However, as in 24], for the energy density 



of the gauge field to increases during the course of inflation we consider / = e c P 2 / nM p with c > 1. 
For our specific example with n = 2 we consider the gauge coupling 



/ 



= cp 2 /2M2 



(75) 



During the first two inflationary stage one can neglect the right hand side of Eq. (|13p and the 
gauge field evolution is given by 

~2 " 



A = p A exp 

where pa is a constant of integration defined in 
in 



-a 



cp_ 

M 2 
p 



(76) 



241 ] . Plugging this into the inflaton equation results 



(p 2 Y + am p 



Acp 2 A 



exp 



.2 i 



-4a 



cp_ 

M 2 
p 



0. 



(77) 



During the first inflationary stage, when the third term in Eq. (|77p is negligible compared to 
the second term, the solution is given by Eq. (|74p . The phase change happens when the last term 
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FIG. 5: Here we plot the evolution of p(a) in chaotic inflation. The dashed-dotted red curve is the analytical 
solution whereas the solid black curve is the full numerical result. As can be seen, the agreement between 
them is very good. The left figure corresponds to e = 0.1 whereas for the right figure e = 0. As explained 
below, the position of the first kink is independent of the value of e and is well approximated by Eq. (|83[) . 
However, the total number of e-foldings depends logarithmically on e. Other parameters are m — 10 Mp, 
p in = 11.2Mp and c = 2.5. 



above becomes comparable to the second term. Because of the attractor mechanism the solution 
rapidly converges to 

/ cp 2 \ m 2 M 2 (c-\) 
ex P " 4 « " 172 • (78) 



M 2 J c 2 p\ 

Note that this solution is analogous of Eq. (|39[) for the symmetry breaking example. Now plugging 
Eq. (USD into Eq. ([77]) yields 

9 

TTt 

3ap= — — p . (79) 

This should be compared to the p equation during the first inflationary phase where now m 2 — > 
m 2 /c. This sudden change of mass induces a jump in p which clearly can be seen in Figj5] . 

Plugging solutions (|78|) and (|74p into Eq. (|76p one finds that A scales like e^ c ^ a and e 3a 
during the first and second inflationary stage respectively. The behavior of the gauge field during 
the first and second inflationary stages has a profile very similar to the plot on the left hand side 
of Fig. [3l It is also instructive to look into the anisotropy parameter. During the first inflationary 
stage one has 

S~|Ri in exp[4(c-lH , (80) 



22 



whereas during the second inflationary stage it reaches the attractor value 

2M 2 c - 1 c - 1 



3p 2 



3c 



(81) 



As expected, S ~ e during the attractor regime. 

Note that Eq. (I77j) can be solved with the general answer 



p 2 ~ 



4Mj,a Mf, 



In 



4(c-l) 



if 



-4(c-l)a 



K - 4(c - l)e c ^/ M p) 



(82) 



where if = Ac 2 p 2 A /m 2 Mp. This has Eqs. (P74|) and (|78l) as the two limiting solutions. In Fig. 
[5] we have compared this analytical solution with the full numerical results and the agreement 
between them is very good. The time of the first phase change, ati, is when the second term in the 
denominator above becomes comparable to the first term which results in 



C P 2 in 



+ 



4(c-l) 4(c-l 



In 



m 2 M 2 (c- 1) 
c 2 p\ 



1 



4( c -i; 



In 



c - 1 M 2 

c 2 



(83) 



This indicates that the smaller is the value of the initial anisotropy Rn n , the longer it takes for 
the system to enter the attractor regime. We have checked that this analytical expression gives a 
good estimate of ol\. 

Like in our previous examples, as the gauge field increases exponentially during the second 
inflationary stage the right hand side of Eq. (|13p becomes important and one enters the final 
inflationary stage. The gauge field equation (in the slow-roll limit) has the same form as Eq. (I44p 
where now the parameter j3 is given by 

6e 2 M P (c- 1) 



c 2 p\ 



(84) 



As in symmetry breaking case the solution is given by the Bessel function Eq. (|49p . The start of the 
third inflationary stage, a = is when the argument of the Bessel function becomes comparable 
to unity so one obtains 02 — -r In (3. Our numerical analysis shows that this expression gives a 
good estimate of 02- As the argument of the Bessel function increases exponentially, the gauge 
field starts to oscillate rapidly. This in turn produces an oscillating effective mass for the inflaton 
in the form of e 2 p 2 A 2 e~ 2a and the slow roll conditions are terminated quickly, ending inflation 
abruptly. Our numerical analysis shows that usually inflation ends when the gauge field makes 
one or two oscillations in less than one e-fold. The behavior of the gauge field during the final 
inflationary stage is similar to the plot on the right hand side of Fig. [3] for the symmetry breaking 
potential. Also the behavior of the inflaton field as a function of time is presented in Fig. [6l The 
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FIG. 6: Here we present the evolution of the inflaton field with respect to the time coordinate in chaotic 
inflation. The left figure corresponds to e = 0.1 whereas for the right figure e = 0. In this plot, the existence 
of the second kink and the duration of the third inflationary phase can be seen clearly when we turn on the 
gauge coupling e. 



start of the third inflationary stage and the existence of the second kink can be seen clearly when 
we turn on e. However, as mentioned above, the third inflationary stage is very short, less than 
an e-fold. This can also be seen from Fig. [5] where the evolution of the inflaton field is presented 
as a function of e- foldings. 



VI. CONCLUSION AND DISCUSSIONS 



In this work we studied anisotropic inflation in models with charge scalar fields. We have 
shown that, similar to 



241 ] . the system reaches the attractor solutions sometime during inflation 
where the ratios 5/e and Ri/e, measuring the level of anisotropies, become at the order of unity. 
This attractor mechanism is fairly independent of the initial conditions. One can tune the model 
parameters such that the time of phase changes, denoted here by a± and ct2, take place within the 
first few e- foldings relevant for the CMB. 

The new interesting effect in our model is the effect of the gauge coupling e on the dynamics 
of the inflaton field and the gauge field. At the final stage of inflation, the term on the right 
hand side of Eq. (|13|) becomes important and gauge field become highly oscillatory. Because of 
the interaction e 2 p 2 A 2 , the oscillations of the gauge field induce an effective time-dependent mass 
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term for the inflaton field and inflation ends shortly after the gauge field starts to oscillate. We 
have studied these effects in examples of symmetry breaking and chaotic inflation models. Both of 
our main results here, that is the existence of the attractor solutions and the oscillatory behavior 
of the gauge field at the end of inflation, show up similarly in these two models. 

We also studied the effect of the gauge field on the dynamics of the waterfall field in charged 
hybrid inflation model. Because of the coupling e 2 p 2 A 2 the onset of waterfall phase transition 
can be significantly different than in standard hybrid inflation. Furthermore, the highly oscillatory 
behavior of the gauge field and its coupling to the inflaton field can play important roles in the 
studies of tachyonic preheating and reheating. Also, as noticed in [3], the inhomogeneous end of 
inflation can have interesting effects for non-Gaussianities in this model. 

In this work we considered only the background dynamics. It would be interesting to study 
the cosmological perturbations in our backgrounds and compare the results, such as the spectral 
index, the mixing between the scalar and tensor modes and non-Gaussianities, with observations. 
We would like to come back to these questions in a future publication. 



Note added: while this work is completed, we became aware of the work [34[ which has some 



overlaps with our results here. We thank the authors in 



341 ] for letting us know about their work. 
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